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[ Abstract. The stable version of the anomaly-induced inflation does not need a fine tun- 

| ing and leads to sufficient expansion of the Universe. The non-stable version (Starobin- 

CN | sky model ) provides the graceful exit to the FRW phase. We indicate the possibility of 

! the inflation which is stable at the beginning and unstable at the end. The effect is due 

pJL^ , to the soft supersymmetry breaking and the decoupling of the massive sparticles at low 

OO . energy. 

m ■ 

> : 

(SI ■ Introduction 

£^ ' The inflation solves many problems of the Early Universe and there are not so much doubts that 

it really took place (see, e.g. On the other hand, the conventional inflaton-based approach 

requires an exact fine-tuning of the form of the inflaton potential or initial data. Some people 
Q H ' believe that there must be a natural mechanism for inflation, which should originate from the 

Q_i| vacuum quantum effects of matter fields. The desired solution would not require a fine-tuning 

f-H | neither for initial conditions nor for the graceful exit to the Friedmann-Robertson- Walker (FRW) 

power-low expansion phase. The purpose of this article is to suggest a qualitative version of such 
a mechanism. Our approach is based on the Starobinsky model [Q, ||, |], ||, ^] and supersymmetry 
(SUSY). We shall also make use of the fact that SUSY is not seen in the low-energy phenomena. 

Consider the vacuum quantum effects in the Early Universe, when the typical energy of quantum 
processes is very high but below the Planck scale. Then, the appropriate framework is not the string 
theory but some quantum field theory. Furthermore, at the energies greater than the masses of 
the particles one can apply an approximation in which the masses of the fields are negligible. The 
matter filling the Universe is characterized by pressure p and energy density p. The standard 
relation p = 3p holds in the ultra-relativistic limit, consequently, the matter decouples from the 
conformal factor of the metric. This means local conformal (Weyl) invariance for the quantum 
fields ||. Then, the vacuum quantum effects of the matter fields can be taken into account through 
the conformal anomaly [||, ||, ||, ||, ||] and may lead to inflation ||. On the other hand, the masses 
of particles should be relevant at lower energies. At the intermediate scales, the anomaly-induced 
action can serve as a model for the mass-independent part of the effective action [||, ^, while 
the leading effect of particle masses would be the renormalization of the Einstein-Hilbert term in 
the vacuum action. Taking this renormalization into account, we arrive at the tempered form of 
inflation g. 

1 On leave from Tomsk Pedagogical University, Russia. E-mail: shapiro@fisica.ufjf.br 
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The paper as organized as follows. In the next section we present the qualitative description of 
the inflation and the graceful exit mechanism. In section 3 some numerical estimates will be done, 
in particular we establish the relation between gravitational scale and the radiation temperature. 
The last section contains brief general discussion of the anomaly-induced approach. 



2. Quantum effects, inflation and graceful exit 

Consider the gauge theory which has No real scalars, Aq/ 2 Dirac spinors, and N\ vectors. We 
shall be mainly interested in the vacuum effects and suppose that the interaction between quantum 
fields is weakened at high energies due to the asymptotic freedom. Than the one-loop vacuum 
contributions play the leading role. Notice that the vacuum quantum effects originate from the 
virtual particles, therefore the numbers iVo,i/2,i do not describe the real matter which might fill 
the Universe. 

The classical action of vacuum can be very complicated and include infinitely many local and 
non-local terms. However, in order to meet a renormalizability requirement there must be, at least, 
the following three terms || [l(]] : 

S vac = J d A x^g~ {hC 2 + l 2 E + l 3 DR } . (1) 

Here, l\ 2 3 are some parameters, C 2 is the square of the Weyl tensor and E is the integrand of 
the Gauss-Bonnet topological invariant. Action (||) consists of the conformal invariant and surface 
terms. Therefore, if we add (|l|) to the Einstein-Hilbert action, then the usual (homogeneous and 
isotropic) FRW solution remains unaltered. 

The renormalization of the action (|l]) leads to the conformal anomaly [11, ||, |l0| 

<T£>= - {wC 2 + bE + cUR) , (2) 

where w, b, c are the /3-functions for the parameters 1%, I2, h 

N + 6iV 1/2 + 12iVi _ _ N + llJV 1/2 + 62iVi _ N + 6N 1/2 - I8JV1 

W ~ 120 ■ (4tt) 2 ' 360 • (4vr) 2 ' ° ~ 180 • (4tt) 2 ' ( } 

One can consider the cosmological model directly, using the anomaly ||, however, it is useful to 
construct the anomaly-induced effective action. The quantum correction V to the classical action 
of vacuum is related to the anomaly 
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The solution of this equation can be found explicitly in the form [12, |T^ [: 

d 4 xV^{waC 2 + ba(E - -OR) + 2ba Aa} ^— J d 4 x^R 2 , (5) 

where we use, along with the original metric g^, another variables: the new metric with the 
fixed determinant and a, such that g^ v = g^ v ■ e 2a . The effective action (||) includes an unknown 
conformal- invariant functional S c [g^ u ], which is the only indefinite component of the solution (||). 
This term is indeed irrelevant when we consider the cosmological (homogeneous and isotropic) 
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metrics. In this case (|5|) is the exact one-loop correction to the classical action of vacuum. The 
second term on the r.h.s. of @ is non-covariant but it can be rewritten in a covariant non-local 

form g g. 

Now we are in a position to discuss the arbitrariness in the classical action of vacuum. It 
is important that this is the action of an external gravitational field and that the metric is not 
quantized. Then, as it was already mentioned above, one can add to the vacuum action ([l]) any 
local or non-local terms. In particular, one can introduce terms similar to those which emerge as 
quantum corrections (||). In this way, one can change the numerical coefficients of all terms in 
anomaly (||) or even cancel the anomaly completely. Of course, if we insist that the classical action 
of vacuum should be local, then only the last J ^—gOR-term in (||) can be modified by introducing 
the / y/— gR 2 -term into the classical action. Sometimes, this operation is called "introducing finite 
counterterm" . 

The introduction of some ( / ^/— gR 2 -type or non-local) extra terms into the classical action 
of gravity may be equivalent to the introduction of the inflaton-like fields (see, e.g., 1€, [l7|]). 
From our point of view, the necessity to adjust the classical action of vacuum after the calculation 
of quantum corrections means that the program of "natural" inflation fails. Hence, we are not 
going to introduce any special vacuum terms. Neither the coefficient of the J ^—gR 2 -teicm in the 
classical action of vacuum will not be taken unnaturally large. On the contrary, we assume that 
this coefficient is essentially smaller than the one in (||) and then, for the sake of simplicity, set it 
to zero. 

The cosmological model is based on the action 

Stotal = — : n ^ [ d 4 X \/—g R + Smatter + S vac + f , (6) 
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where the quantum correction V is given by (^). It proves useful to introduce the following 
variables: conformal factor a = e a , physical time t (where dt = a{rj)dr) and rj is the conformal 
time) and H{t) = &(t). For the conformally flat case (similar solutions for the FRW metric with 
k = ±1 are also possible, see |J) g^ u = n^, the equation for H(t) has the form Q 

H + 7HH + A (3-^ HH 2 + \H - ± b -H A - — L_ ^2H 2 + h) = , (7) 

There are special solutions with H = const: 

H = and H 1/2 = dr — ==== = , a(t) = a ■ exp Ht . (8) 

y — lovrGo 

Indeed the positive sign in the last expression corresponds to inflation Q. 

The next step is to study the stability of inflationary solution under fluctuations of the conformal 
factor a{t). The analysis of stability can be performed analytically Q or numerically j6|. Solution 
(|8|) is stable under fluctuations of H (t) if the particle content of the quantum theory satisfies the 
condition b/c < [[§. By using Eq. (||), we arrive at the following criterion of stability [17]: 

Ni < i N 1/2 + -L N . (9) 



2 The non-local nature of the anomaly contribution has been first noticed in If £ 
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One can imagine that there was some string phase transition at the Planck scale. When all massive 
string modes decouple, we are left with the ultra-relativistic matter described by field theory. Then, 
if @ is satisfied, the inflation starts independent on the initial conditions. The numerical analysis 
of H shows that the stabilization of inflation performs in a fraction of the Planck time and the 
Universe needs just about 100 Planck times (we define the Planck time as tp = 1/Mp = G l l 2 ) 
to expand into necessary 65 e- folds. The main problem is that the stable inflation is eternal and 
no simple receipt is known for the graceful exit to the FRW phase. The possible solutions of this 
problem were discussed in Refs. || [l^] using the effective field theory approach. In particular, 
we have supposed that when the typical energy decreases, the masses of matter particles become 
relevant and the deviation of a(t) from the exponential behavior might lead to the graceful exit. 
To some extent, this letter is devoted to the realization of this idea (see also M). 

Let us also mention Ref. [O]. In this interesting paper the classical action of vacuum is chosen 
in such a way that the total / y/— gR 2 -teim in (^) is absent. Then the inflation is stable but it can 
be destabilized by metric fluctuations. Indeed, this kind of solution is theoretically possible, but 
we will not discuss it here. 

The non-stable version of the anomaly-induced inflation g, [7|, ||, [HJ has another advantage: it 
easily breaks and one can achieve the graceful exit to the FRW phase M. Furthermore, there can 
be rapid oscillations of the conformal factor after the breaking, potentially producing the reheating 
H). The effective action (g) includes the massive degree of freedom associated to a. The decay of 
this mode into matter particles has been discussed in |3[ [/], [|, 19]. We remark that this massive 
mode is not a fundamental field but instead it is a degree of freedom induced by the quantum 
effects of matter fields. If some matter field decouples (say, because it has a large mass), then it 
does not contribute to the massive mode of a anymore and the coefficients Nq, Ny2, N% change 
correspondingly. At this moment the energy of the induced mode can transfer into the real (not 
virtual) matter sector through the creation of particles. In the time scale, the last decoupling is 
the one of the lightest neutrino. After that the particle content in (||) is N\ = 1 (one photon) 
and Nq = N^m = 0. It is easy to see that in this case the inflationary solution is unstable while 
the solution Hq in (||) is stable. Due to this circumstance, in the present day Universe we do not 
have fast inflation! This simple example shows, by the way, the predictive power of the anomaly- 
induced inflation. The existence of more than 18 massless scalars or more than 6 Weyl fermions is 
completely ruled out by the condition (|9|). 

The shortcoming of the non-stable inflation is that without some " strong measures" it does not 
last long enough and the Universe does not expand sufficiently. The necessary "strong measures" 
have been discussed in S, [5| . They consist in the extremely exact fine tuning of the initial conditions 
and in the introduction of the / yJ—gB? term into the classical action ^. 

Let us summarize. We have two sorts of the anomaly-induced inflation: stable and unstable. 
The advantage of the stable one is that it does not depend on the initial data. However, it remains 
unclear how the inflation stops. On the contrary, the non-stable inflation stops immediately and one 
is forced to use the "strong measures" in order to achieve the necessary expansion of the Universe. 
Indeed, both versions do not look perfect if they are considered separately. The best situation 

3 I am very grateful to A. A. Starobinsky for explaining to me this point. 
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would be to have an inflation which is stable at the beginning and unstable at the end. If one could 
switch from one to another in a natural way, this could be the desired explanation of inflation. 

Now we come back to the condition of stability (^), which depends on the particle content. 
The known spectrum of particles fits with the Minimal Standard Model (MSM). Let us present a 
few details. The SM includes 6 quarks, each of them has 2 chiralities and 3 colors. Hence, quarks 
contribute N quar k s = 18. Furthermore, there are 6 leptons. Taking all the neutrino massive, we 
arrive at Ni eptons = 6. The scalar sector has one Higgs doublet Ao = 4, while the vector one 
consists of 8 gluons, W^, Z and photon. After all, we have Ao = 4, Aq/2 = 24, Aq = 12 . Then, 
Eq. @ indicates to the non-stable inflation and this perfectly fits with our dream to have unstable 
inflation at the end. The same concerns, obviously, the present-day Aq = 1, Ao = Aq/ 2 = case 
considered above. 

Now, let us remember that the anomaly-induced inflation is supposed to occur at the sub-Planck 
energy domain. There are many reasons to expect that the particle spectrum at this scale goes 
beyond the SM. In particular, it may happen that the high-energy theory possesses super symmetry. 
Let us, for example, look at the particle content of the Minimal Supersymmetric Standard Model 
(MSSM). This model has Aq = 12 as in the SM. Besides the known fermions, one needs additional 
superpartners (gaugino) to all the vectors. The same concerns Higgs particles, which require 
higgsino. In total, we have Aq/ 2 = 32 for the MSSM. Finally, the scalar sector includes two Higgs 
doublets and numerous superpartners of the fermions: squarks and sleptons, such that A^o = 104 . 
It is easy to see that this particle content provides stability in (^). In fact, similar result can 
be expected for any realistic supersymmetric model. The supersymmetric extension of the gauge 
theory implies the replacement of any vector multiplet by the N = 1 vector superfield. Moreover, 
one has to add superpartners to the fermions in such a way that they form chiral superfields. Both 
operations increase the number of spinors and scalars in (||) and we arrive at the stable inflation. 





Figure. The samples of the one-loop Feynman diagrams which contribute to the anomaly. 
The bubble of the matter field has 2,3 or 4 external lines of the field a. 

Thus, the situation when the anomaly-induced inflation is stable at the beginning and unstable 
at the end means exactly that supersymmetry breaks at the last stage of inflation. If we suppose 
that the typical energy scale decreases during the inflation, this can be associated with the super- 
symmetry breaking at low energies. Let us explain the last statement in more detail. The Feynman 
diagrams which contribute to the anomaly-induced action (|5|) consist of a quantum bubble of mat- 
ter (non-gravitational) fields with external tails of the a field (see the Figure). According to the 



Appelquist and Carazzone theorem [20|, the loop of massive field decouples when the energy of 
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external lines becomes much smaller than the mass of the quantum field in the loop. One has to 
notice that if the origin of masses of quantum fields is not the Spontaneous Symmetry Breaking, 
the decoupling theorem applies nicely. The massive sparticles decouple when the typical energy 
of the external lines of the field a becomes smaller than the masses of these particles. Therefore, 
the graceful exit is realized most naturally for the soft SUSY breaking which is, also, the most 
acceptable from the phenomenological point of view (see, e.g. [2l]]). The mechanism of decreasing 
the typical energy of gravitons during inflation will be discussed in the next section. 



3. Some numerical estimates 

It is important to estimate the duration of inflation until the SUSY breaks. Since the breaking 
of the unstable inflation occurs in a very short time ||, this is equivalent to the two questions: 

i) What is the SUSY breaking scale? 

ii ) How to evaluate the energy of the cr-field quanta in the external lines of the diagrams? 
The first problem has been widely discussed in the literature (see, e.g. pi] ). The most popular 

is the situation when the supersymmetry breaks softly above the electroweak scale Mp « 300 GeV . 
There are no upper bounds for SUSY, and its breaking may occur at the GUT scale Mx ~ 10 16 GeV 
or even at the Planck scale M P = G~ l l 2 = 10 19 GeV, depending on the SUSY model. 

To answer the ii) question is not a simple problem. One can try to use the framework of 
renormalization group in curved space-time |jTo| . This version of the renormalization group links 
the scaling of all dimensional quantities with the one of the metric. Then, as any dimensional 
quantity, the typical energy of the cr-quanta should vary as fj, a ~ l/o, while the metric transforms 
as 9fiu 9av • a 2 - Unfortunately, this consideration is consistent only for the constant scaling 
parameter a. In the case of Eq. (|8|) the dynamics of such important dimensional quantity as scalar 
curvature is different from the above rule, and in fact it is constant R = — 12H 2 . Since H and R 
are the most important local dimensional parameters associated to a(i), we can identify the energy 
of the a-quanta with H = a /a. At low energies, when the high derivative terms in (|5|) become 
negligible, this also follows from the Einstein equations which we shall apply below. 

After identifying the graviton energy with H we meet another problem. The decoupling mecha- 
nism described above works if the energy scale is decreasing during inflation, but in the exponential 
phase (||) the Hubble parameter H(t) = H\ does not decrease. However, this is true only if we do 
not take the masses of the quantum fields into account. Compared to the massless case, the fourth 
derivative anomaly-induced terms are still present in the effective action of massive fields. Besides, 
there are three other contributions. First of all, one meets the renormalization of the Newton 
constant G, and its scale dependence. In a parallel paper we have shown that this effect can be 
taken into account and results in the following evolution of H: 



where iVj and rrii are multiplicity and mass of the fermion of the specie i. 

Another important effect is the renormalization of the cosmological constant. In this paper we 
do not consider this effect which will be discussed in details elsewhere. It is sufficient to mention 
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that the cosmological constant does not change the graceful exit mechanism. Finally, there are 
contributions given by the infinite power series in curvature. Unfortunately, there is no available 
method for calculating these terms in a general form. Some existing approximations may be 
misleading, in particular because we can not use the weak curvature limit. Let us suppose that 
in the high energy region the leading effect is the renormalization of G and perform our analysis 



using Eq. (|i0|) 



For simplicity we suppose that the SUSY breaks at GUT scale, and that the gauge group and 
particle content are such that 

£ 3(^)2 N i m i = (WMxf = W 34 GeV 2 = l(T 4 Mj, . (11) 

Since the numbers iVo, Ni/2,Ni are much greater than for the MSSM, the magnitude of the Hubble 
parameter Hi will be essentially smaller than ]-[( MSSM ) ~ Mp. The exact relation between Hi 
and Mp requires fixing the particle content of the SUSY GUT. Without going into these details, 
let us suppose that in this situation Hi = Mp/10. The inflation will last until the Hi(a) becomes 
comparable to Mx- Then ( |l0|) and ( |il~l) immediately lead to the relation 

a(Hi = M x ) « a • e 100 . (12) 

It is easy to see that our suppositions about the spectrum of GUT fermions were not in favor of 
too much expansion, but still we arrived to much more than the minimal necessary 65 e-folds. The 
lower energy scale for the SUSY breaking can increase the number of e-folds to enormous extent. 

The next important problem is to link the energy of the a(i)-quanta with the temperature T r 
of the electromagnetic radiation. Consider the lower energy region after the SUSY breaks down 
and the inflation becomes unstable. Then the higher derivative terms in the effective action are 
irrelevant and the graviton energy H can be defined from the Einstein equation with the radiation 
dominating the matter stress tensor. The standard estimate ]l]] shows that 

T r w (H ■ Mp) 2 . (13) 

Taking H ~ Mp we obtain T r ~ Mp that corresponds to the unification of all fundamental sources 
and the limit of applicability of the semiclassical approximation. Taking (as in our estimate leading 
to ([l2])) the value H = Mx, we meet the temperature satisfying Mx < T r < Mp. In order to 
illustrate the gravitational suppression of T r , let us consider the energy when the lightest neutrino 
decouples from gravitons H = m u ~ 10~ 3 eV. According to (13) the corresponding temperature 



is T r > 10 3 GeV . This temperature is much higher than the QCD and EW scales, such that all 

vectors of the MSM should be considered massless. According to (Jq|) , the instability of inflation is 

guaranteed. An important consequence is that the nucleosynthesis can not be jeopardized by the 

anomaly-induced inflation. 

Another aspect of the energy scale problem is the following. The change of the regime of 

expansion due to the decoupling of massive sparticles should provoke yet another quantum effects 

4 This is nice, because our semiclassical approximation (no quantum gravity, string etc) becomes really consistent. 
At lower energies Hi in could be greater, but this does not matter at all, because after the SUSY breaking the 
exponential solution gets unstable and inflation does not occur. 
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like particle creation. The new particles absorb the energy of the decaying massive mode of the a 
field, as it was discussed in []|, ||. Also, during the inflation (||), the real matter filling the Universe 
is out of the equilibrium state and does not lose energy. Therefore, the content of the Universe after 
the end of inflation are both hot matter which came from the string epoch and did not completely 
cool down, and the hot matter created during the transitional period from the stable to the unstable 
inflation. After this transition is over, the Universe entered into the FRW phase and all processes 
proceed in a standard way. 

After the transition into an unstable phase, there can be, in general, very different versions 
of further behaviour for a(t) (2|, [2^, ||. The analysis of the phase structure |J shows that a(t) 
has various attractors and only some of them can be associated to the graceful exit to the FRW 
solution. The evolution of a(t) phase can be studied numerically. The result depends on the initial 
data a(0), d(0), d(0), d(0), which are not well defined since we do not have the detailed description 
of the decoupling. In this situation it is better to consider different versions of these initial data. 
For example, if we take the purely massless case (0) and substitute the initial data a^ k \0) = H\ 
with Hi of the MSSM) into the equation with the MSM parameters b and c, the system goes to 
the undesirable "hyperinflation" ||. However, if we take the masses into account, the situation 
changes drastically. From the numerical analysis of || follows that at the end of the stable period 
the evolution goes in a power-like manner a(t) ~ t l l n . Then the initial data for the second MSM 
phase must be taken such that each derivative must be much smaller than the previous one. The 
typical time of the stable phase is 100 tp. If we take, using the Planck time units f| 

o(0) = 1 , a(0) = 10" 2 , a(0) = 10" 4 , 'd'(O) = 10~ 6 , (14) 

then the asymptotic behaviour of a(t) is quite similar to the desirable a(t) = i 1 / 2 . The possible 
small changes in ( |T4] ) do not modify the qualitative behaviour of a(t). Of course, these results must 
be accepted with the proper caution. Let us remind that by the end of the decoupling the approxi- 



mation (10) is not safe and one has to take into account other contributions to the effective action. 
This can lead to the more precise definition of the initial data a' 4 '(0) for the unstable inflation. 
The subsequent progress in the anomaly-induced inflation requires quantitative description of the 
decoupling process (unfortunately, this is a difficult problem similar to the one people meet in QCD) 
and of the transition to the FRW phase. After that the metric and density perturbation analysis 
(see, e.g., |7|, ^3], |l^, should provide useful information for the phenomenological investigations. 



4. Discussions 

It is interesting to discuss the general status of the anomaly-induced inflation. Theoretically, 
the phenomena of anomalous breaking of the local conformal symmetry is well established (see, 
e.g |1H|). The same concerns the existence of the inflationary solution (|8|). Of course, there is no 
guarantee that this kind of inflation really took place in the Early Universe. However, there is only 
one way to avoid the anomaly-induced inflation: to introduce the classical action of vacuum in such 
a way that the inflation becomes unstable at all scales. As an example we mention the y/—gR 2 -term 
which could provide inequality b/c > independent on the particle content Nq, Nim,N\. Then, 
the inflation must have another origin. But, in this paper we have shown that the absence of extra 
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vacuum terms may produce the natural inflation and the natural graceful exit due to the SUSY 
breaking. At this level, there are no contradictions neither with the present-day Universe nor with 
the nucleosynthesis mechanism. 

Further investigations of the density and metric perturbations are necessary. They can brink us 
to the point when experiments and observations can show whether the anomaly induced inflation 
really took place. 
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